VT symmetry with a system of three-level atoms 
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We show that a vapor of multilevel atoms driven by far-off resonant laser beams, with possibility 
of interference of two Raman resonances, is highly efficient for creating parity-time (VT) symmetric 
profiles of the probe-field refractive index, whose real part is symmetric and imaginary part is anti- 
symmetric in space. The spatial modulation of the susceptibility is achieved by proper combination 
of standing-wave strong control fields and of Stark shifts induced by a far-off-resonance laser field. As 
particular examples we explore a mixture of isotopes of Rubidium atoms and design a PT-symmetric 
lattice and a parabolic refractive index with a linear imaginary part. 
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While non-Hermitian operators obeying pure real spec- 
tra, like for example the Bogoliubov-de Gennes [l[ spec- 
tral problem, linear stability problem for nonlinear waves, 
or simply a parabolic potential with linear imaginary 
part [2J, are known in physics for long time, it was only 
due to the work Q that fundamental importance of such 
operators became widely recognized. It was discovered 
in Q that there exists a wide class of complex poten- 
tials of Schrodinger equation obeying pure real spectra, 
and even most importantly, that this property is intrinsi- 
cally related to the parity (V) and time (T) symmetries 
of physical systems. This discovery triggered the discus- 
sion [4j on the fundamentals of quantum mechanics whose 
axioms are based on Hermitian operators for observables. 
Further growth of interest in the theory of parity-time 
(VT) symmetric potentials was originated by suggestions 
of implementation of VT symmetry in a waveguide with 
gain and absorption [f| , which was based on the analogy 
between quantum mechanics and paraxial optics where 
the refractive index plays the role of the potential in the 
Schrodinger equation. In optics PT-symmetric refractive 
index profiles (i.e. obeying gain and losses of a given ge- 
ometry) has been experimentally realized using process 
of four- wave mixing in Fe-doped LiNb03 substrate [6J. 
The possibility of optical realization of VT symmetric 
potentials motivated various suggestions of practical ap- 
plications, like nonreciprocal wave propagation 0,0] , im- 
plementation of coherent perfect absorber [§| , giant wave 
amplification Q, etc. Experimental realization of PT- 
symmetry using plasmonics [lOj and temporal simulation 
of lattices using optical couplers ll| were also reported. 

In this Letter we demonstrate the possibility of 
practical implementation of spatially distributed 'PT- 
symmetric refractive index, i.e. the one having the prop- 
erty n(x) = n*(—x), in vapors of multi- level atoms driven 
by control fields with properly chosen Raman resonances 
and by a far-off-resonant Stark field. 

First we recall some recent achievements in creation of 
large susceptibilities in atomic vapors controlled by ex- 



ternal laser beams. While such systems are intrinsically 
dissipative, it was suggested in [12J and shown experi- 
mentally in (l3| . that using the destructive interference 
in imaginary part of susceptibility, it is possible to achieve 
large real refractive indexes keeping the absorption small 
enough. Recently, the idea of using two far-off-resonant 
control fields for realizing high susceptibility with nearly 
zero absorption of a probe field was developed theoreti- 
cally 14 [ and confirmed experimentally (TBj . An alterna- 
tive way of achieving similar effect was proposed in 16 1 
where two A-systems were explored in order to excite two 
Raman resonances (see also Ref. [17|). 

Since the above mentioned schemes use the interfer- 
ence of two Raman resonances, one of which resulting 
in gain and another in absorption, the imaginary part 
of probe-field susceptibility appears as a non-monotonic 
function of the frequency with both positive (gain) and 
negative (absorbing) domains. Even more remarkable 
property is the possibility to design distributions where 
real and imaginary parts of the susceptibility appear re- 
spectively as even and odd functions of the probe-field 
frequency 14Ul7|]. Since for a monochromatic beam the 
change uj — > —lj is equivalent to the change t — » —t the 
last property can be viewed as a time inversion symme- 
try. Thus our goal can be formulated as completing this 
time symmetry by the symmetry in the coordinate space. 

To achieve this goal we choose the scheme based on 
mixture of two species of A-atoms, similar to the one ex- 
plored in [l6| . The involved atomic states will be assigned 
as \g, s) (ground state), \a, s) (lower state), and |e, s) (ex- 
cited state), hereafter s = 1, 2 indicates the specie of the 
atoms (Fig. [T]). The species have atomic densities Ni 
and N 2 , respectively (N = Nx + N 2 being total atomic 
density) . Vt p is the half Rabi frequency of the probe field 
coupling \g, s) f-> |e, s) and Qi^ are half Rabi frequencies 
of two control fields coupling \a,s) -H> \e,s). All fields 
are far-off resonant, which is guaranteed by condition 
A s 3> S) where A s = cj e>s — w a . s — Wc is one-photon 
detuning with Hu!i, s (I = g, a, e) being the eigen energy 
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FIG. 1: (color online) (a) Two A-systems and Raman transi- 
tions used for obtaining 'PT-symmetric refractive index. Ini- 
tially populated levels are indicated by the filled circles, (b) 
Possible geometry for the suggested scheme. All notations are 
defined in the text. 



of the state \l, s), and lo p (ui c ) being the center frequency 
of the probe (coupling) field. Notice that one can ex- 
plore different cases by changing signs of the two-photon 
detuning, defined by S s = u> a ,s — ^g.s — (w p — w c ). In 
Fig. []Ja) we show an example of the proposed scheme 
where Si > and $2 < 0, what corresponds to the VT 
symmetric lattice obtained below. In our configuration, 
the first scheme (s = 1) exhibits two-photon absorption 
for the probe field while the second one (s = 2) provides 
two-photon gain. 

Spatial modulation of the probe-field susceptibility can 
be achieved by using x-dependent control-field [with Rabi 
frequency Q s (x)]. Such filed, however affects both one- 
and two-photon detunings. Therefore we explore the sec- 
ond possibility, which is modulation of relative energy- 
level shifts along x-direction resulting in the dependence 
A s = A s (x). This task can be achieved if a strong, far- 
detuned laser field Es(x) cos(wst), where E$ and uis are 
respectively amplitude and frequency, is applied to the 
system. This field originates Stark shifts of levels by 
AEs, s (x) — —jai. s E^(x), where a^ s is the scalar polar- 
izability of the level \l,s). Below Es(x) is referred to as 
the Stark field. If within the required accuracy one can 
consider a 3lS sa o. a ,si i-e. the difference of Stark shifts 
between the ground-state sublevels is negligible (this is 
the situation we consider below). Then 5 S is not affected 
by the Stark field, while for the one-photon detunings we 
have A s (x) = A s - (a e , s - a S)S )I5g(x)/(47i). Thus the 
Stark field allows one to manipulate the spatial distribu- 
tion of A s (x), on the one hand, and on the other hand 
being far-off-resonant does not lead to power broadening 
of lines due to low rate of transitions [18J . 

Since the Stark shift usually appears in the second 
order of the perturbation theory, it requires relatively 



strong electric fields. Taking into account the charac- 
teristics of the available lasers, this may impose strong 
limitations, because the characteristic scale of the spatial 
modulation of A s (x) is of order of the wavelength of the 
Stark field As- For the typical order of the control field 
amplitudes E c ~ 10 2 V/cm (fi 1)2 ~ 2?r x 514 MHz), the 
required amplitude of the Stark field must be three or- 
ders of magnitude larger, i.e. E$ ~ 10 5 V/cm (for more 
detail estimates see below). Being focused into a spot of 
a diameter of 30 /im this requires laser powers of order 
of 100 W. Nowadays such powers can be achieved using 
say quantum cascade lasers 19|, 2(J operating at micron 
wavelengths (in [19] it was As = 4.45 /im, which will be 
used in our estimates). 

For the described model the susceptibility of the probe 
field can be computed from the density-matrix formalism 
within the rotating-wave approximation. Its functional 
form is the same as in [l6| with the difference that now 
the two-photon detuning and control field depend on the 
spatial coordinate: 
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Here xo = Nid 2 g 1 /(eoh), Eq is the vacuum permittivity, 
r\ = N%d\ g 1 characterizing the ratio between the 

specie densities is considered as a free parameter, the d.. iS 
stands for the dipole moment of the transition between 
the ground an excited states of the s-the system. 

Notice that generally speaking in a warm vapor 
large Doppler broadened line widths (typically Alo/lu ~ 
10~ 6 [18]) may degrade the effectiveness of resonant 
schemes. While in our case such broadening would not 
be important for large one-photon detuning it is of or- 
der of the two-photon detuning. It turns out, however, 
that the effect of this line broadening can be significantly 
suppressed (by the factor \u c — wJ/w c ) if far-off-resonant 
co-propagating beams (see e.g. [2lJ) with the frequencies 
close enough are used. This is the case we consider below. 

Now our task is to determine spatial distributions of 
f2 s (x) and A s (x) ensuring the condition 



f(x) = n(x) — n*(x) = 0, 



(2) 



where n(x) = n r (x) + irii(x) ss [1 + ^Xp( x )]V^o * s ^ ne 
refractive index, for x being either arbitrary or belonging 
to some interval (see below). u(x) will be used to control 
the accuracy of the obtained 'PT-symmetric refractive 
indexes. 

Because of the complexity of the relations among all 
parameters involved in Eqs. (JTJ, ([2]), it is not obvious a 
priori that the problem has a solution. Therefore we pro- 
ceed with analysis of particular systems. To this end we 
adopt the approach as follows. First, we define a "seed" 
shape of the susceptibility, of the type we would like to 
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obtain, say X \ x i ej, rf), which contains a number of free 
parameters ej and r\. Second, neglecting all terms within 
some accuracy, say of 10%, from © we compute ana- 
lytical solutions for OJ (x, ej, rf) and A sd (x, ej, rf). Third 
we substitute the obtained "seed" Rabi frequencies and 
Stark shifts in Eq. ([1}. Due to crude approximations 
made the so obtained susceptibility may still give signifi- 
cant errors in (j2J . Our final step is to minimize the error 
using the control parameters ej and rj. 

First, we apply the above algorithm to obtain a VT- 
symmetric lattice. Since the main limitation for the 
lattice period is determined by the Stark field, we re- 
quire the period to be As, i.e. x( x ) — x( x + ^s)- 
The experimental geometry we bear in mind is illus- 
trated in Fig. [TJb). We further assume that each con- 
trol field consists two almost parallel plane waves having 
cc-components of the wavevector equal to ks = 2w/Xs, 
i.e. to the Stark field wavevector. Now the seed solution 
can be chosen as x sd = e o + £iCos£ + 162 sin £, where 
£ = ksx and £0,1,2 are free parameters. As a particular 
atomic vapor we use a mixture of isotopes 87 Rb (s = 1) 
and 85 Rb (s = 2), and assign \g,s) = \5S 1 / 2 ,F = 1), 
\a,s) = \5S 1/2 ,F = 2), and \e,s) = \5P 1/2 ,F = 1) for 
each specie. 

To perform numerical study we further particularize 
the problem by assuming that Ni = 1.26 x 10 12 cm -3 
and N 2 = 1.14 x 10 12 cm" 3 (thus N = 2.4 x 10 12 cm" 3 ). 
The isotopes are loaded in a cell at approximately 363 K 
and the coherence decay rates are estimated as ^ eg pa 
7 eQ = 2?r x 334 MHz and ~/ ag = 2tt x 16 kHz [lej. With 
sufficiently high accuracy we can impose a e ,i — ~ 
a e ,2 — a g ,2 — 27rfix 0.1223 Hz(cm/V) 2 and d ae ,i ~ d ae ^ — 
2.5377 x 10~ 27 C-cm 22], what allows us to consider 51i ps 
fl 2 = Sl c - Other parameters are chosen as lo p ~ lo c = 
2ttx3.77x10 14 s" 1 {\ p = 795 nm), rj = 0.91, xo = 0.57 s, 
61 = 1.807e S , and 5 2 = -0.0l7 es . 

As the second step, we solve equation x P (£l| d , A s s d ) = 
X sd with respect the real ttf and Af. For |<S 2 | < \Si\ 
it is possible to leave only the leading terms in S 2 al- 
lowing one to find the solutions explicitly. As the final 
step, we substitute £l s s d and A s s d back into Eq. (fTJ) and 
obtain v(x) as small as possible by tuning the parame- 
ters £0,1.2- Since the obtained formulas have rather cum- 
bersome forms they are not presented here. Instead, we 
write down the final expressions for the control-field Rabi 
frequency and for the Stark-field amplitude by keeping 
the first significant harmonics (i.e the terms > 10~ 4 ) 



■ £ [0.9698- 
7e 9 [1.5384- 



0.0053 cos £- 
-0.0122 cos£ 



0.0007 sin £], 
-0.0232 sin }]. 



Here Eq = 4 x 10 5 V/cm [see Fig. |5{a)]. As we already 
mentioned, the final shape of the susceptibility differs 
from the ansatz x sd - We represent it also in a form of a 
Fourier series by keeping the first significant harmonics, 
i.e. Xp( x ) — Xop+Xipi x )i where xop ~ 0.2257 determines 
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FIG. 2: (a) The Rabi frequency of the control (solid line) 
and Stark (dashed line) fields vs. £. (b) n r and (c) rii vs. £. 
(d) Real (solid line) and imaginary (dashed line) parts of the 
error function 



the average (in space) refractive index and the periodic 
part is given by 

Xi P « 0.0075 cos^ + il0" 4 [3.9418sin^ + 0.2298 sin(2£)]. 

. (3) 

The real and imaginary parts of the refractive index is re- 
spectively shown in Fig. [2jb), (c). We notice that there 
is a significant difference in the magnitude of the con- 
stant real part of the susceptibility and its imaginary 
part. What is important, however, is that the latter con- 
stitutes about several percents of the variation of the real 
part. In Fig. [2fd) we show errors in obtaining the real 
and imaginary parts of the susceptibility. 

Now turn to the propagation of the probe field along 
z-dircction in the described atomic vapor. Since in the 
y-direction the medium is homogeneous, in the paraxial 
approximation the beam propagation is governed by 



2i/up „ 
oz 



d 2 n p 

dx 2 



klxip( x )tt p = 0, 



(4) 



where k p — to p ->Jl + Xop/c is the probe- field wavevec- 
tor. [The beam of constant amplitude in y-direction is 
chosen only for convenience: one can consider any waveg- 
uide structure in y-direction what will result only in the 
renormalization of the constants in Q]. By the ansatz 
fl p (x, z) = £l(x)e lbz where b is the propagation constant, 
we obtain 



B = ^-b 



¥ h - (5) 



Taking into account the well known results on the peri- 
odic potentials 23] as well as rapid decay of the Fourier 
harmonics in ([3]), one expects that with high accuracy 
the spectrum of ([5]) is indeed pure real. To check this, 
in Fig. 0Ha),(b) we show Q obtained numerically for the 
susceptibility including all terms until the leading term 
violating the condition @ (the latter have amplitudes 
< 10- 5 ). 
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FIG. 3: (Color online) Real (a) and imaginary (b) parts of the 
dimensionless propagation constant j3. Solid (dashed) lines 
show the results for PT-symmetric refractive index and non 
PT-symmetric refractive index obtained by change of the mu- 
tual concentration of species, as explained in the text. Evo- 
lution of |f2 p (2;)| 2 in the vapors with PT-symmetric (c) and 
non "PT-symmetric (d) refractive indexes. 
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FIG. 4: Real (a) and imaginary (b) parts of the refractive 
index. Solid line is the refractive index inside the atomic cell 
limited by Bragg mirrors. Dashed line shows the distribution 
if the gas would occupy the whole space, (c) Spatial distri- 
bution of Stark (dashed line) and control (solid line) fields 
required for inducing the local PT-symmetric refractive in- 
dex, (d) Real (solid line) and imaginary (dashed line) parts 
of the error function inside the cell. 



Let us now clarify whether the described procedure is 
"structurally stable", i.e. whether small deviations of 
system parameters do not break the obtained PT sym- 
metry. To this end we tested the method with respect 
to change of mutual concentration of the species. More 
specifically we have changed the obtained n by 10% and 
repeated the calculation described by the above algo- 
rithm (i.e. we did not use 77 as a matching parameter 
any more) . We indeed found that now the accuracy with 
which @ is satisfied is lower (see Fig. [3fa),(b) ) but the 
imaginary parts are still very small (of order of 10~ 3 



In Fig. [3] we also illustrate propagation for the input 
Gaussian beam in the discussed 'PT-symmetric [panel 
(c)] and non- 'PT-symmetric [panel (d)] structures. The 
input probe beam £l p (z = 0) = e -0 - 1 ^ 3 *) propagates in 
the PT-symmetric medium along much longer distance 
compared to the non PT-symmetric case where absorp- 
tion is observed. 

The described approach can be used to produce other 
shapes of the refractive index. As an example we now 
show how the method can be modified to obtain a PT 
symmetric parabolic refractive index Q ■ The main idea 
is based on the fact that PT-symmetry can be satisfied 
only locally in space. Then one can "cut" undesirable 
non-PT-symmetric distribution of the refractive index 
by choosing the location of the finite-size vapor cell with 
respect to the domain where the external field produces 
local 'PT-symmetry. This is illustrated in Fig. H] (a),(b) 
where we still consider the mixture of the Rubidium iso- 
topes, but taking now 85 Rb as the first, s = 1, and 87 Rb 
as the second, s = 2, systems. We choose Ni — 2.39 x 
10 15 cm" 3 and N 2 = 8.98 x 10 12 cm" 3 , rj = 0.375 x 10~ 2 , 
Si = 3.56 x 10~ 4 7 eff , and 5 2 = 3.93 x 10 _3 7 eg , without 
changing the other parameters. The described procedure 



of the choice of the Stark and control fields is performed 
to satisfy ([2]) only locally in space. In particular, this 
is achieved, by taking O c = j eg (2 + 0.2£ + 0.009£ 2 ) and 
E s = T o (0.3695 - 0.13750, illustrated in Fig.^c). The 
spatial modulation of the susceptibility now has a com- 
plex form [Fig.@]a, b]. If however the atomic cell (limited 
say by Bragg mirrors) is situated as shown in the pan- 
els, then the atomic vapor changes the refractive index 
only in the domain occupied by the cell. This domain 
was chosen to ensure @ with high accuracy where the 
real [solid line in panel (a)] and imaginary [solid line in 
panel (b)] are respectively parabolic and linear. Now 
the susceptibility inside the atomic cell is described by 
X P = 0.2021 - 0.0007c; 2 - i0.0006£, which obvioulsy sat- 
isfies the condition of PT symmetry. 

To conclude, we have suggested a possibility of creat- 
ing PT-symmetric profiles of the refractive index in a 
mixture of resonant atomic gases. An important prop- 
erty of the proposed scheme is the possibility of chang- 
ing parameters of the structure in situ, say by changing 
wavelengths or amplitudes of the control fields. While we 
used the well studied isotopes of the Rubidium atoms the 
proposed scheme allows for further generalization and im- 
provement using other atomic isotopes, four-level atoms, 
mixtures of more than two isotopes, monoatomic vapors 
with two control fields, etc. Adding more control param- 
eters opens possibilities of experimental implementation 
of nonlinear PT-symmetric susceptibilities 25j, as well 
as combined linear and nonlinear ones [26| which recently 
attracted much attention. 
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